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I. IHTRODUCTIOK 


The interaction of ocean surface waves with a permeable, 
rubble mound breakwater is complex. In general, the 
incident waves are partially transmitted through or over the 
porous structure, partially dissipated by breaking on the 
rough slopes and by turbulent friction within the 
breakwater, and partially reflected. The reflection of 
waves from breakwaters (and other reflectors such as sand 
bars and seawalls) is usually measured with an array of 
pressure sensors or surface height gauges deployed seaward 
of the reflector (e.g., Thornton and Calhoun 1972; Mansard 
and Funke 1980; Yokoki et al. 1992). The interpretation of 
these measurements is complicated owing to the fact that the 
incident and reflected waves are phase-coupled. 

In many studies the array analysis is simplified by 
assuming that the incident waves are uni-directional, 
propagating perpendicular to the reflector. Thornton and 
Calhoun (1972), Morden et al. (1976) and Goda and Suzuki 
(1976) used two sensors positioned on a line perpendicular 
to the reflector to decompose the wave field into incident 
and reflected wave contributions. This technique breaks 
down at the frequency where the wavelength is equal to twice 
the sensor spacing. Mansard and Funke (1980) overcame this 
problem by applying a least-squares-fit technique to three 



sensors, and this approach was extended to linear arrays 
with an arbitrary number of sensors by Zelt and Skjelbreia 
(1992) . 

The assumption of normally incident waves used to 
estimate reflection is often violated in a natural coastal 
environment. Refraction of swell over complex bathymetry 
may result in significantly oblique angles of incidence at 
the breakwater, especially if the breakwater is not aligned 
with the surrounding depth contours. Directionally broad, 
locally generated seas are typically only weakly refracted 
in depths greater than 10 m and can approach a breakwater at 
relatively large oblique angles. Although in principal 
array measurements seaward of a reflector can be used to 
infer the reflection of a directionally spread wave field, a 
very large number of sensors is required to obtain reliable 
estimates of the directional spectra of both incident and 
reflected waves (Isobe and Kondo 1984) . 

In this study, a new method is presented for estimating 
wave reflections from array data. The assumption of normal 
incidence used in previous studies is relaxed, but angles of 
incidence are assumed to be small. Reflection from the 
breakwater is assumed to be a linear process governed by 
Snell's law (i.e., specular reflection). An expansion for 
small oblique incidence angles is used to derive approximate 
relationships between the array cross-spectra and a small 
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number of parameters that describe the incident wave 
properties and the reflectivity of the breakwater. 

The new estimating technique is applied to pressure 
array data acquired offshore of the Monterey Harbor 
breakwater at Monterey, California. The field site, array 
geometry and data acquisition are described in Section 2. 
Section 3 presents the new estimation technique, and 
demonstrates the accuracy and robustness of the method 
through model simulations. Results of the analysis of 
reflections observed at the Monterey breakwater are given in 
Section 4, followed by a discussion and conclusions in 
Section 5. 



II. EXPERIMENT 


The experiment was conducted at the Monterey Harbor 
breakwater in Monterey, California (Fig. 1). Situated at 
the southern end of Monterey Bay (Fig. la), the permeable 
breakwater is exposed only to swell with deep water 
propagation directions ranging from west to north-northwest 
and seas generated by local northerly winds. The 
(refracted) directional wave spectra at this sheltered site 
are typically narrow with most of the energy approaching the 
breakwater from northerly directions within about 20° of 
normal incidence to the breakwater. 

Construction of the 400 m long rubble mound breakwater 
was completed in 1934, and that same year a 120 m long 
extension was added (Fig. lb). During subsequent years, 
much of the interior of the breakwater washed away and the 
elevation of the structure decreased so that it was often 
overtopped at high tide (Thornton and Calhoun 1972) . In 
1993, the breakwater extension was improved by increasing 
the elevation back to the original design height of +3 m 
relative to MLLW, and by adding on the harbor side a single 
cover layer of 8-ton-average armor stone. A cross section 
of the breakwater can be found in Thornton and Calhoun 
(1972) . The present field study was conducted at the 
improved extension. 



A 2-dimensional array of six pressure sensors was 
deployed seaward of the breakwater in a depth of about 16 m 
(Fig. lb) to measure both the incident and reflected wave 
fields. The small depth variations (less than 1.5 m) across 
the array are neglected. The array aperture is 31 m x 15.6 
m and the distance between the breakwater and the closest 
sensor P6 is about 52 m. Additionally, a single pressure 
sensor P7 was positioned inside the harbor at a distance of 
about 47 m from the breakwater to measure the transmitted 
wave field. Sensor locations relatively close to the 
breakwater were chosen to avoid contamination by edge 
effects (e.g., diffraction around the tip of the 
breakwater). 

The array data were acquired with a tattletail micro¬ 
processor located in the center of the array, which was 
cabled to a radiotelemetry system on the breakwater, from 
where the data were transmitted to a computer at the Naval 
Postgraduate School. Data were acquired with a 2 Hz sample 
rate continuously for eight months between August 1993 and 
March 1994. 

Typical power spectra measured along a line 
perpendicular to the breakwater, converted to surface height 
with a linear theory depth correction, are shown in Fig. 2. 


The swell spectrum is most energetic in the frequency range 
of about 0.06-0.1 Hz. Large differences in energy levels 



is the vector wavenumber of the reflected wave. R, (|> and k,. 
are generally functions of k and O). Adding Eqs. 1 and 2 
yields the total surface elevation function T|(x,t): 

Ti {X, t) = Jjexp [ i (k-x-w t) ] «3 > 

+R exp[i{ic^' x-(i)t+(l)) ]] dz{(ji,k) 

The cross-spectrum h^((0) of two sensors at positions x;, 
and 3c„ is defined as: 

J exp [iwt] i3n„(u) du = £{t| {x„, t) t|* (x„, t+t)} (4) 

where E{} denotes the expected value and * the complex 
conjugate. Substitution of Eg. 3 in Eg. 4 yields: 

+R exp [i {k^- x^+^) ] ] 

X [exp[-ikxj +R exp[-i(JCj,' x„+4)) ] ]5(«,k) d* 
where S{( 0 ,k) is the wavenumber-freguency spectrum: 

5(o) k) - (<«>,*) dz* (<0 < Jic)} (5) 

dcodk 

The wavenumber magnitudes of both incident and reflected 
waves are given by the linear dispersion relation 
ctf=gk[tanh(kh) ] , where k=|k |=|kr|, g is gravity and h the 
water depth. Introducing polar coordinates 
(k=[k(cos(0)),k(sin(0))]) and further assuming that the 




is specular (kr=[-k(cos(0)),k(sin(0))]), Eg. 5 



can be expressed as (dropping the frequency dependence): 

J^G^(0)5{0) d0 (7) 

where G^^O) is given by: 

G^{ 0 ) =[exp [iic(x„cos 0 +y„sin 0 ) ] 

+i?{ 0 ) exp [i (-icx„cos 0 +icy^in 0 +(j){ 0 )) ]] 

X [exp[-i^:(x„cos 0 +ynSin 0 ) ] 

+i?( 0 ) exp [-i (-icx„cos 0 +icy„sin 0 + 4 ) ( 0 ) ) ] ] 

and S{0) is the directional spectrum of incident waves. 

Assuming that incident waves approach the breakwater at 
small oblique angles, R and (|) can be expanded for small 0: 

i?{0) =i?(,+i?20^+i?40^+. . . (9) 

(t){ 0 ) =<l)o+<t)20"+<t)40''+. . . (10) 

The odd-order terms in these expansions vanish because R(0) 

and 0(0) are symmetric functions. Substituting Eqs. 9 and 
10 in Eq. 8 and expanding G,„(0) for small 0 yields : 

with the lowest three coefficients A^, and 0 ,^^ given by: 

A^= exp [iic(x^-x,,) ] + 2 i?gCOs [k(x„+x„) - 4 ) 5 ] 

+i?ifexp [-ik(x„-x„) ] 
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B^=ik(y^-y^)A^ 


(12b) 


^Vm-yr) )exp [ ik {x^-xj ] 


+ [ 2 i? 2 -j?oic 2 {y„-y„) 2] cos{k{x„+xJ 
+ [ 2 J?„<j) 2 +i?gic(x„+x„) ] sin(A:{x^+x„) -<|)(j) 

+ [|i?o'(-ic2 (y^-y„) 2 +iic(x„-x„))+ 2 i?gi? 2 ] exp [-iic(x„-x^) ] 


(12c) 


Ga^id) can be accurately approximated by a truncated 
expansion of the form Eq. 11 if all the expansion 
coefficients (Eq. 12) are <0(1). Thus the technique 
presented here is generally valid only for compact arrays 
(i.e., |k • x|<0(l)) positioned close to the reflector 
(i.e., within a wavelength so that <|>S0(1)). 

The exact G„„,(0) (Eq. 8) of sensor pair n=3, m=l (Fig. 
lb) for f=-;^=.08 Hz (the dominant swell frequency) is 
compared to the linear (A^+Bn„6) and quadratic (A^+Bnj„0+C„j,0^) 
approximations in Fig. 3. In this calculation the 
breakwater is idealized as a partially absorbing wall at 
x=L=83 m with R independent of 0 and <l)=2kL [cos (0) ] (Fig. 
lb) . The linear approximation diverges from the exact 
solution for |0|>5-15° but errors are generally within 10% 
even for incidence angles as large as 30°. The quadratic 
approximation is more accurate for small values of 0, but 
diverges sharply from the exact solution for large 0 with 
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deviations exceeding 20% for |0|>3O°. Similar calculations 
for other sensor pairs and other frequencies in the swell 
range (0.05-0.12 Hz) show comparable agreement between the 
exact G^(0) and the linear and quadratic approximations. 

Substitution of the truncated expansion of G^(0) (Eqs. 
11,12) in Eq. 7 yields the quadratic approximation of the 
cross-spectrum h^^: 

,4>o ) E+B^ (i?o »‘l>0) 

, , (13) 

+ C^(i?0,i?2,(l)o,4)2) [0Ls+0Lan]^ 

where E, 0„ean 0n„s are the energy spectral density 

(integrated over all directions), mean propagation direction 
and root-mean-square directional spread of the incident 


E=jis{d)de 

f^0S(0)d0 

Q _-'-2 (14 b ) 

fls(e)de 

(d-e^^j^s(d)dd ^ 

•'-I _ (14c) 

fls(e)dd 
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Eq. 13 relates the array cross-spectra to three 
parameters of the incident wave field (E, 0„ean- Qrms) and 
four reflection parameters (Ro , 00 , , (t >2 ) . In general, 

inverse algorithms can be developed that search for a 
combination of these seven parameters on a frequency by 
frequency band basis such that the associated cross-spectra 
are as close as possible to the observed cross-spectra, but 
this requires an extensive array. Since the six-element 
array deployed in the present study is relatively small in 
aperture, the contributions of the quadratic (C^^) terms to 
Eq. 13 is small even for wave incidence angles as large as 
20° (e.g., compare the linear and quadratic approximations 
of GjiO) in Fig. 3). Thus, the higher-order parameters 
Rj, and 02 i^ot be resolvable within the uncertainty of 

the array cross-spectra. Neglecting the quadratic terms in 
Eq. 13 reduces the inverse problem to only four unknowns: 
two incident wave parameters (E, 6„ean) two reflection 

parameters (Rq , 0o) • Estimates 0„ean' $o were 

obtained by minimizing a simple root-mean-square misfit norm 




(15) 


with h^ the observed cross-spectra and h^ the linear model 
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< -^0 ' $o > < -^0 - $0 ) ^mean^ 


(16) 


Since both the number of sensors and the number of unknowns 
is small, this minimum misfit can be evaluated by 
essentially sweeping through the entire parameter space. 

For all possible combinations of E, $o' the 

cross-spectra h,^ and the misfit £ were calculated with Eqs. 
12a, b, 15 and 16 to obtain a global minimum for e. The 
range of physically plausible values for E, 6„ean- snd 
used in these calculations is 0.25E-4E, -30°-30°, 0-1, and 
0-360°, respectively, with 1 the average of the auto spectra 
h„„. The minimum value of £ was calculated by sweeping 
through all possible combinations of E, 0„ean, Rq, <t)o< 
stepping with increments of 0.15E, 5°, 0.1, and 15°, that 
were somewhat coarse owing to limited computing resources. 
The accuracy of the solution fi, 6„eaa, Rq, $o was improved by 
sweeping with smaller step sizes (0.03E, 2°, 0.02, 3°) 
through a reduced parameter range (fi±0.3E, , fto±0.2, 

^o±30°) . 

The accuracy of the new estimation technique was 
verified through a series of model tests with simulated 
array cross-spectra. For a chosen incident swell 
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directional spectrum of the form: 


S(0) 


the true cross-spectra h^ of the Monterey array (Fig. lb) 
were evaluated with Egs. 7 and 8. As before, the breakwater 
was idealized in these tests as a partially absorbing wall 
at x=L=83 m with R independent of 9 and (t)=2kL [cos (0) ] . In 
each test, five random realizations of cross-spectra h^j, 
were generated using the procedure described in Long and 
Hasselmann (1979). Errors in the simulated finm include both 
statistical uncertainty resulting from finite length data 
records (160 degrees of freedom) and uncorrelated instrument 
noise (noise to signal ratio 0.05) and are roughly 
representative of the actual array measurements presented in 
Section 4. The simulated were then treated in exactly 
the same way as cross-spectra obtained from ocean 
observations, and the minimum misfit e and the optimal 
parameters fi, , Ro, were obtained with the inverse 

algorithm described above. 

Results of model tests for typical swell with frequency 
f=0.086 Hz are presented in Figs. 4 and 5. Similar 
comparisons (not shown) for other frequencies in the swell 
band (0.05-0.12 Hz) yielded comparable agreement. Estimated 
(R) and true (R) reflection coefficients are in good 
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agreement for both normally (Fig. 4a) and obliquely (Fig. 

4b) incident swell (discrepancies generally within ±0.05). 
Estimates of the incident wave energy E (Fig. 4c) and 
propagation direction 0„ean (Fig- 5) also agree well with the 
input spectrum (errors generally less than 20% and 4°, 
respectively). The values of the misfit £, a normalized rms 
measure of the discrepancies between the model {hr.j„) and 
"observed" (hn..„) cross-spectra (Eq. 15), range from 0.02 to 
0.08 (Fig. 4d). The misfits are generally larger for the 
simulations with obliquely propagating waves than for 
normally incident waves (e.g., compare the e for 0j„ean=O° 

20° in Fig. 4d) owing to errors in based on a small 0 
approximation (Eq. 16). However, aside from a slight bias 
in R (Fig. 4b), the estimates of incident and reflected wave 
parameters do not appear to be significantly degraded for 
0mean=2O° (FigS. 4,5). 

Overall, the model simulations demonstrate that the 
estimation technique is relatively insensitive to errors in 
the data and the model, and can extract accurate estimates 
of wave reflections from the array data acquired in the 
present study. 
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IV. OBSERVED REFLECTIONS 


The reflection estimation technique described in the 
previous section was applied to the array measurements 
acquired at the Monterey Harbor breakwater (Fig. 1). Nine 
days were selected for analysis, that approximately span the 
range of conditions encountered during the experiment. For 
each of these nine days, two three hour data runs were 
processed, one acquired at low tide and one at high tide, 
with the objective to examine the sea level dependence of 
breakwater reflections. Cross-spectra with a frequency 
resolution of 0.0078 Hz and 160 degrees of freedom were 
computed for each of the 18 data runs. For every frequency 
band in the dominant swell range, 0.05 Hz-0.12 Hz, the 
inverse algorithm was applied to the cross-spectra to 
estimate the incident wave spectral density E(f), the mean 
incident wave propagation direction 6r„ean(f)< the reflection 
coefficient Ro(f) and the phase lag between incident and 
reflected waves (|)o(f). Frequencies less than 0.05 Hz and 
greater than 0.12 Hz are not considered here because 
infragravity waves dominate the spectra below 0.05 Hz (e.g., 
Fig. 2; Okihiro et al. 1992) and above 0.12 Hz (i.e., 
frequencies greater than about twice the spectral peak 
frequency) local nonlinear effects may be significant. 


16 



The variability of incident wave conditions was small. 
Estimates of the incident swell variance ranged from 3.7 cm^ 
to 51 cm^ (i.e., significant wave heights of 7-29 cm). 

These low energy conditions are typical for this site owing 
to sheltering effects (Fig. 1), and the fact that no major 
storm occurred during the eight month data acquisition 
period. Estimates of the mean frequency and propagation 
direction of incident waves (averages of f and 0„,,n(f) over 
the swell band, weighted by E(f)) ranged from 0.065-0.089 Hz 
and from 0-18°. These incidence angles are well within the 
range for which the present technique (based on a small 6 
expansion) is expected to be accurate (Figs. 3-5). 

Estimates of Ro(f) obtained from different data runs are 
remarkably similar (Fig. 6), even though the incident wave 
spectral levels varied by more than an order of magnitude. 
Reflection of small amplitude swell from the breakwater is 
apparently insensitive to the incident wave amplitude, 
consistent with the assumption that reflection is a linear 
process. In all cases the reflection estimates show a 
strong frequency dependence with Ro(f) decreasing 
approximately linearly with increasing frequency from about 
0.7-0.8 for f=0.05 Hz to 0.2-0.3 for f=0.12 Hz. In contrast 
to estimates reported by Thornton and Calhoun (1972), the 
present observations do not suggest a strong dependence of R 
on tidal sea level variations. 



The misfit e between and (Eq. 15) is shown in 
Fig. 7 as a function of frequency. In the most energetic 
part of the spectrum (0.06-0.1 Hz), e is approximately 0.03- 
0.07, comparable to the misfits obtained in model tests. At 
frequencies below 0.06 Hz and above 0.1 Hz, where energy 
levels are relatively low, the misfits are slightly larger 
(0.05-0.16) than expected from model tests, possibly owing 
to directional spreading and/or nonlinear effects. 

Estimates of the phase lag (t»o(f) between incident and 
reflected waves obtained from 18 different data runs are 
approximately equal, increasing with increasing frequency as 
expected from theory (Fig. 8). Neglecting depth variations 
(i.e., changes in k) seaward of the breakwater, the 
theoretical phase lag ()>o for small incidence angles is equal 
to 2kL. This crude approximation of <t)o, taking L to be the 
distance to the crest of the breakwater, is in good 
agreement with the estimates across the entire swell band 
(Fig. 8). 


18 



V. DISCUSSION AND CONCLUSIONS 

A single sensor {P7) was deployed inside the harbor 
(Fig. lb) to obtain crude estimates of wave transmission 
through the breakwater. Assuming that the propagation 
directions of incident, reflected and transmitted waves are 
nearly perpendicular to the breakwater, the fluxes of 
incident (Fi(f)), reflected (Fj.(f)) and transmitted 
energy are approximately given by: 

F,{f)=E{f)Cg,{f) {17a) 

F,{f)=R^{f)E{f)Cg^{f) (17b) 

F,{f)=E^{f)Cg,{f) (17c) 

where Cg^ and Cg^ are the group velocities at the offshore 
array and sensor P7, E is the estimated incident wave 
spectrum and E„ is the transmitted wave spectrum measured 
by sensor P7. The transmission coefficient T(f), defined 
as: 

(18) 

is plotted as a function of frequency in Fig. 9 for all 18 
data runs. The observed transmission coefficients are more 
variable than the reflection coefficients (Fig. 6), but do 
not suggest a consistent sea level or frequency dependence. 
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Estimates of the residual energy flux : 

Fa{f)=F^(f)-F^(f)-F^(f) ( 19 ) 

that is dissipated through wave breaking on the rough 
breakwater slope and/or turbulent friction inside the porous 
structure, are variable, but F,j/Fi generally increases with 
increasing frequency (Fig. 10) . 

Estimates of the bulk incident, reflected, transmitted 
and dissipated energy fluxes, obtained by integrating Fi(f), 
Fj.(f), Ft(f) and F,j(f) over the entire swell band (0.05-0.12 
Hz), are presented in Fig. 11. The observed fraction of the 
incident energy flux that is reflected from the breakwater 
(Fig. 11a) varies between about .2 and .5. These changes in 
the breakwater reflectivity are primarily the result of 
variations in the dominant swell frequency (0.06-0.09 Hz, 
Fig. 6). On the other hand transmission of wave energy 
through the breakwater appears to be a strong function of 
the incident energy flux (Fig. 11b). On days with very low 
amplitude swell about 40-60% of the incident energy flux is 
transmitted through the breakwater and dissipation is weak 
(0-40%, Fig. 11c). On more energetic days about 40-60% of 
the incident energy flux is dissipated (Fig. 11c) and the 
transmission is reduced to about 20-30% of the incident 
energy flux (Fig. 11b). 

Estimates of energy transmission and dissipation also 
show a dependence on sea level. At high tide the 
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transmitted (dissipated) energy fluxes are slightly larger 
(smaller) than at low tide, possibly owing to the fact that 
the effective width of the breakwater near the sea surface 
is smaller at high tide than at low tide. 

Waves transmitted through the breakwater may undergo 
partial reflection at the other side of the harbor and the 
associated standing wave patterns may contribute significant 
errors to estimates of F, (f) based on Eg. 17c. Accurate 
estimation of wave transmission requires an array of sensors 
on the harbor side of the breakwater, which was not 
available in this study. Furthermore, the range of 
conditions encountered in the present experiment was rather 
limited, and the reflection/transmission coefficients 
observed under benign conditions may not be representative 
for the breakwater performance under storm or large 
amplitude swell conditions. More extensive measurements are 
needed to evaluate the performance of permeable rubble mound 
breakwaters. 

The main result of the present study is the development 
of a new technique for estimating the reflection of a 
random, directionally spread wave field from a coastal 
structure (e.g., a breakwater or seawall) or natural sand 
bars. The estimation technique can be applied to a compact 
(i.e., aperture less than a wavelength) array of pressure 
sensors or surface height gauges of arbitrary geometry. 


21 



deployed seaward of the reflecting surface (within a 
wavelength). Model tests demonstrate that for wave 
incidence angles less than about 30° the new method can 
provide accurate estimates of the gross properties of 
incident and reflected waves. 
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APPENDIX 


Figure 1. a) Location of the Monterey Harbor breakwater, 
b) Coordinate frame and locations of pressure sensors 
deployed near the 120 m breakwater extension. The three 
sensors in the center (P2,P4,P5) form an equilateral 
triangle with dimension 2.5 m. PI is 15 m seaward of the 
triangle, P3 and P6 are both positioned at a distance of 13 
m from the center. Sensor P7 is inside the harbor. 
Soundings are in meters relative to MLLW. 

Figure 2. Typical power spectra of array sensors PI 
(solid), P5 (dashed), P6 (dotted), positioned on a line 
perpendicular to the breakwater (Fig. lb). The dash-dot 
curve is the spectrum measured inside the harbor (sensor 
P7) . 

Figure 3. The exact G,^(6) (solid lines, Eq. 8) compared to 
quadratic (left panels) and linear (right panels) 
approximations (dashed lines, Eq. 11), for sensor pair n=3, 
m=l, with the reflection coefficient, R, equal to 0, .5 and 

1 . 

Figure 4. Results of model tests for normally (0^30=0°) and 
obliquely (0maati=2O°) incident waves, with the reflection 
coefficient, Ro, varying from 0 to 1. a) Rq vs R^ for 
0mean=O°- t*) Rq VS Rq for Q^^^=20°. c) Ratio between 
estimated (E) and true (E) incident wave energy vs R,,. d) 
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between the estimates and true values. Pluses and circles 
in panels c) and d) correspond to 9 „„.„ =0° and 20° 
simulations. 

Figure 5. Results of a model test for a reflection 
coefficient Ro=. 5 and mean incidence angles 0„ean varying from 
0°-20°. The estimated incidence angle is compared to 

9mean- The solid line denotes perfect agreement. 

Figure 6. Estimated reflection coefficient Ro vs. frequency 
for all 18 data runs. The solid and dashed lines indicate 
low and high tide runs. 

Figure 7. Misfit e between h^ and h^ as a function of 
frequency for all 18 data runs. The solid and dashed lines 
indicate low and high tide runs. 

Figure 8. The estimated phase lag between incident and 
reflected waves (relative to x=0) as a function of frequency 
for all 18 data runs. Circles and crosses indicate low and 
high tide runs. The solid line is a crude theoretical 
approximation of ([»(,. 

Figure 9. The transmission coefficient T as a function of 
frequency for all 18 data runs. The solid and dashed lines 
indicate low and high tide runs. 

Figure 10. The ratio of dissipated to incident energy flux 
as a function of frequency for all 18 data runs. The solid 
and dashed lines indicate low and high tide runs. 

Figure 11. The ratios of reflected to incident (a), 
transmitted to incident (b), and dissipated to incident (c) 
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Figure 11. The ratios of reflected to incident (a), 
transmitted to incident (b), and dissipated to incident (c) 
energy fluxes vs. the incident energy flux for all 18 data 
runs. The energy fluxes are integrated across the swell 
band (0.05-0.12). The circles and crosses represent low and 
high tide runs. 
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